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ABSTRACT 

Aims. We investigate the nature of the classical low-velocity structures in the local velocity field, i.e. the Pleiades, Hyades and Sirius moving 
groups. After using a wavelet transform to locate them in velocity space, we study their relation with the open clusters kinematically associated 
with them. 

Methods. By directly comparing the location of moving group stars in parallax space to the isochrones of the embedded clusters, we check 
whether, within the observational errors on the parallax, all moving group stars could originate from the on-going evaporation of the associated 
cluster. 

Results. We conclude that, in each moving group, the fraction of stars making up the velocity-space overdensity superimposed on the 
background is higher than the fraction of stars compatible with the isochrone of the associated cluster. These observations thus favour a 
dynamical (resonant) origin for the Pleiades, Hyades and Sirius moving groups. 
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1. Introduction 

It has been known for a long time that the local velocity field 
in the solar neighbourhood is clumpy, and that most of the ob- 
served clumps are made of spatially unbound groups of stars, 
called moving groups (e.g., Eggen 1958, 1960, 1975, 1983). 

Although unbound, most of these moving groups how- 
ever share the kinematics of well-known open clusters. The 
best documented low-velocity groups (e.g., Famaey et al. 
2005; Chereul et al. 1998, 1999; Dehnen 1998; Montes et al. 
2001; Ecuvillon et al. 2007) are the Hyades moving group 
(e.g., Famaey et al. 2007) associated with the Hyades cluster 
(600 Myr) and the Sirius moving group associated with the 
evaporating UMa star cluster (300 Myr). Another kinematic 
group called the Local Association or Pleiades moving group 
is a set of stars mostly associated with the Pleiades cluster 
(100 Myr), and with a few other young clusters with ages all 
below 150 Myr (a Persei, NGC 2516, IC 2602 and Scorpio- 
Centaurus). Finally, the large group consisting of stars lagging 
behind the galactic rotation and moving outward in the disk 
(e.g., Blaauw 1970, Raboud et al. 1998, Famaey et al. 2005, 
Bensby et al. 2007) is called the Hercules stream. 
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The classical hypothesis ofEggen is that those low-velocity 
moving groups are in fact remnants of clusters which partly 
evaporated with time. Indeed, their prominence close to the 
Sun's position in velocity space without any vertical flow (e.g., 
Seabroke et al. 2007) is not compatible with the hypothesis that 
they are the result of merger events with satellite galaxies. Such 
satellite debris have been identified in the solar neigbourhood 
(Helmi et al. 2006, Arifyanto & Fuchs 2006, Dettbarn et al. 
2007), but containing much less stars and moving with higher 
velocities. However, Eggen's hypothesis is still debatable be- 
cause at least some of the moving groups may also be generated 
by resonant mechanisms linked with the non-axisymmetry and 
non-stationarity of the Galaxy. E.g., a rotating bar at the cen- 
tre of the Milky Way is now thought to be able to create the 
Hercules stream if the Sun is located near the outer Lindblad 
resonance (OLR) of the bar. The existence of such a kinemat- 
ical group could be due to the coexistence at the OLR of or- 
bits elongated along and perpendicular to the bar's major axis 
(Dehnen 2000), or to an overcrowding of chaotic regions in ve- 
locity space induced by the bar near the OLR (Fux 2001). On 
the other hand, moving groups with lower velocities such as the 
Hyades and the Pleiades could still be linked with the spirality 
of the Galaxy. E.g., a series of strong transient spiral arms with 
their mean corotation at the solar galactocentric radius have 
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been shown to produce similar structures in the local velocity 
distribution (De Simone et al. 2004), while a group such as the 
Hyades could also correspond to nearly-closed orbits trapped 
at the 4 : 1 inner Lindblad resonance (ILR) of a two-armed 
spiral density wave (Quillen & Minchev 2005). 

Eggen's hypothesis and the dynamical perturbations by the 
bar and the spiral arms are of course not mutually incompatible 
scenarii. Clusters do evaporate over time and, at some inter- 
mediate stage, an unbound group of stars with similar veloci- 
ties must appear (Woolley 1961). On the other hand, large spi- 
ral perturbations are known to indeed exist in the Galaxy, and 
must have an effect on local stellar kinematics. So, for an indi- 
vidual moving group, how can we observationally distinguish 
between a pure Eggen scenario, a pure dynamical scenario and 
a mix of both? 

Famaey et al. (2005) used a maximum-likelihood method 
(Luri et al. 1996) to model a sample of K and M giant stars 
(combining Hipparcos parallaxes with Tycho-2 proper motions 
and CORAVEL radial velocities) in order to assign each star 
to a kinematic base group. This allowed to identify three main 
moving groups (Hercules, Sirius and Hyades-Pleiades) and to 
plot the corresponding stars in Hertzsprung-Russell (HR) dia- 
grams: the similarity between the HR diagrams of stars in the 
moving groups and in the background was noted, which ar- 
gued in favour of a dynamical (resonant) origin for the groups 
rather than a coeval (evaporated cluster) origin. However, this 
method, although extremely powerful in order to yield unbi- 
ased distance estimates, presented two shortcomings: first, the 
method revealed itself unable to separate distinct substructures 
on a scale smaller than 30 km/s (e.g., to separate the Hyades 
from the Pleiades moving group); second, although the sim- 
ilarity of HR diagrams for the moving groups and the back- 
ground was striking, it did not statistically exclude the possi- 
bility of such HR diagrams being the result of coeval groups 
superimposed on the background. However, a complementary 
analysis of the Hyades moving group, by Famaey et al. (2007), 
based on the Geneva-Copenhagen catalogue of F and G dwarfs 
(Nordstrom et al. 2004), revealed that the number of stars in 
the group with respect to the total number of stars in the back- 
ground was independent of mass, contrary to the expectations 
based on the difference between the initial mass function (IMF) 
of an evaporating cluster and the present-day mass function 
(PDMF) of the background, arguing again in favour of a res- 
onant origin for the Hyades moving group. 

Here, we reanalyze the data of Famaey et al. (2005) in 
order to circumvent the two shortcomings of the maximum- 
likelihood method. First, we use a wavelet transform (see 
Chereul et al. 1998, 1999 who had not used the radial veloci- 
ties that are now available) in order to identify the structures on 
scales smaller than 30 km/s. Second, by directly comparing the 
location of moving group stars in parallax space to isochrones 
of the embedded clusters, we independently confirm the dy- 
namical nature of the groups (Famaey et al. 2005, 2007) in a 
robust statistical way. 



2. Data 

Fifteen years of observations with the CORAVEL spectrome- 
ter on the Swiss 1-m telescope at the Observatoire de Haute 
Provence yielded the radial velocities of Hipparcos stars later 
than about F5 in the northern hemisphere. This unique database 
of radial velocities measured with a typical precision of 
0.3 km/s, combines a high precision and the absence of kine- 
matic bias. The details of the sample of K and M giants used 
in this paper (combining Hipparcos parallaxes with Tycho-2 
proper motions and CORAVEL radial velocities) can be found 
in Famaey et al. (2005). 

3. Locating the moving groups 

The parametric maximum-likelihood technique of Luri et al. 
(1996) and Famaey et al. (2005) is not suited to find small 
structures in the velocity space and cannot separate for example 
the Hyades and Pleiades moving groups. To detect these small- 
scale structures, we use a wavelet transform technique on the 
2-D velocity field defined by the U and V velocities (U is the 
velocity towards the galactic center, V the velocity in the di- 
rection of Galactic rotation, both with respect to the Sun). The 
method is similar to Chereul et al. (1998, 1999) analysis but for 
clarity, the main steps are summarized below. 

Let us however note that the method used in Famaey et al. 
(2005) had the great advantage of giving an unbiased statistical 
estimate of the distance of individual stars, unlike the inverse 
of the parallax affected by the Lutz-Kelker bias. In the pro- 
cess of analyzing the distribution of stars in velocity space with 
the wavelet transform, we thus have to make use of the results 
of Famaey et al. (2005) concerning the distances of individual 
stars. 

The wavelet transform by an oscillating and zero-integral 
mother wavelet *P(x) of a real one-dimensional signal F (x) is 
defined as: 

W s (i) = T J F(x) ■ V* (-^-jdx (1) 

where s is the scale and i the position of the analysis. The 
values W s (i) are called the wavelet coefficients. The original 
signal F(x) can then recovered by a double summation of 
W s (i) • ¥(— ) over the variables s and ;. 

Following Chereul et al. (1998, 1999), we use the a trou al- 
gorithm (e.g., Holschneider et al. 1989) to compute the wavelet 
decomposition of our 2-D velocity field. This algorithm relies 
on the definition of the mother wavelet v F(x) as the difference at 
two different scales of the same smoothing function (or scaling 
function) O(x): 

Y(x) = 0(x)-i<D(|), (2) 

and on the relation 

^(^^Z^W^-O. (3) 

Here h(l) is a one-dimentional discrete low pass filter ap- 
plied iteratively to compute the smoothed distributions C s at 
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scale s. As in Chereul et al. (1998, 1999), we use h(l) = 
^hhbh^ 0Tl = (-2,-1,0, 1,2), and <D(jc) = Qx - 2| 3 - 
4|x - 1| 3 + 6\x\ 3 - 4\x + 1| 3 + \x + 2| 3 )/12. 

For our 2-D analysis, we define <b(x,y) = <f>(x) ■ 0(y). The 
coefficients C s (i, j) [i.e. the smoothed signal at scale s and pixel 
(i, j)] are given by the recursive relation 

1=2 m=2 

C s (j, J) = Yj X ^ )/l(m) C - l(/ + 2 " 1/ ' j + 2 "' m) ' (4) 

/=-2 m=-2 

Co being the observed distribution ,F(t/, V) convolved by the 
scaling function 0(jf,y) (see Chereul et al 1999, Fig. 1 for a de- 
tailed schematic overview of the procedure). With this relation, 
the distance between two bins increases by a factor 2 between 
the scales s and s+ 1. The wavelet coefficients W s (i, j), are then 
given by the signal difference between two scales 

W s (i,j) = C s - l (i,J)-C s (i,j). (5) 

It is clear from this relation that the wavelet coefficients 
at each scale have the same spatial coverage in the velocity 
space as the original distribution, which permit to localize 
the structures in the velocity space. In practice, we bin the 
2-D velocity field F(U, V) on a 500x500 squared grid from 
-250 km/s to +250 km/s. We then analyze the observed 
histogram on 6 dyadic scales s corresponding to 5, 7, 10, 14, 
20 and 28 km/s with the previous definition. 

The original distribution is subject to Poisson noise and 
noise fluctuations induce non-zero wavelet coefficients at each 
scale. These non-zero coefficients generate artifacts mimicking 
substructures and must be removed prior to the scientific anal- 
ysis. To reject the noise-induced coefficients and select only 
the regions corresponding to real overdensities, we again fol- 
low Chereul et al. (1999). We perform numerical simulations 
of a 2-D Poisson noise and follow rigourously the same proce- 
dure (binning and wavelet transform) to estimate the amplitude 
of the wavelet coefficients generated by the noise at each scale. 
We then compute a treshold A* (j, j) at each scale and each pixel 
such that the probability that the noise wavelet coefficient is 
larger than this threshold is P — 10~ 4 . We subsequently apply 
a hard thresholding to clean the wavelet coefficients using 



W s (i,j) 



\W s (i,j) if W s (i,j)>%(i,j) 
otherwise 



(6) 



The small-scale structure appears then very clearly on the 
14 km/s scale, which is plotted on Fig. 1. 

The wavelet coefficients at scale 4 (see Fig. 1) are non- 
zero in the three boxes of the [U, V] -plane corresponding to 
the Pleiades (17 a [-23,-9], V [-28,-18]), the Hyades 
(U * [-41,-23], V * [-24,-12]), and Sirius (U * [0,10], 
V - [-5, 8]). Interestingly the Hercules stream is also re- 
covered, but is separated in two different substructures (U - 
[-65, -49] and [-40 - 30], V [-55, -47]), a fact which had 
not been identified so far in previous studies. 

The most interesting result from this wavelet analysis is 
thus the similarity between the small-scale structure of the 




Fig. 1. Contours and wavelet coefficients in the [U, V] plane 
for the Famaey et al. (2005) sample. Black: isocontours of 
the velocity distribution smoothed over 14 km/s. The color 
scales with the wavelet coefficients at scale 4 (14 km/s) af- 
ter thresholding. The four structures observed in the wavelet 
space correspond to the Hercules stream (bottom green struc- 
tures), Pleiades and Hyades (middle red structures from right to 
left) and Sirius (yellow, top). The limit of the structures in the 
wavelet space are used to select the samples for further analy- 



Hyades-Pleiades moving group and the Hercules stream. They 
both appear as a pair of groups. If they are both of dynamical 
origin, this could provide an important new constraint on dy- 
namical models designed to reproduce them (see e.g. figure 9 
ofDe Simone et al. 2004). 



4. Age of the moving groups 

Now that the locations of the substructures have been precisely 
identified in velocity space, we look whether these subtsruc- 
tures result or not from the sole evaporation of the main clus- 
ters associated with them. Since the dynamical nature of the 
Hercules stream is well-proven from a detailed analysis of its 
chemical heterogeneity (Bensby et al. 2007), we concentrate 
hereafter on the Pleiades, Hyades and Sirius moving groups. 

For each zone in velocity-space determined hereabove with 
the wavelet technique, we shall first estimate the fraction of 
stars in the moving group and in the background. Then, we 
shall assume that all the stars of the group are coeval and thus 
located on a single isochrone in the Hertzsprung-Russell di- 
agram. We can then compute for each star the parallax that 
would correspond to this isochrone. Then, the relative differ- 
ence (normalized by the error) of the Hipparcos parallaxes and 
isochrone ones in parallax space should follow a gaussian of 
mean and standard dispersion 1 after the distribution of back- 
ground stars in parallax space has been subtracted. If this is 
not the case, the moving group will be proven to be, at least 
partially, of dynamical (resonant) origin. 
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4.1. Estimating the fraction of stars in the moving 
groups 

In order to separate the background from the Pleiades compo- 
nent, the histogram of stars in the strip -28 < V (km s _1 ) < -18 
(corresponding to the range spanned by the Pleiades in Fig.Q]) 
was plotted in Fig. |2j and three separate gaussians fitted to 
the total histogram (see Famaey et al. 2007): one represent- 
ing the background, another representing the Pleiades and the 
last one partially representing the Hyades (which partially con- 
taminates the {/-distribution in this strip). Since the Hyades 
V-strip is slightly different from the Pleiades one, the opera- 
tion has been repeated for the Hyades with the strip -24 < 
V (km s~') < -12. Finally, for Sirius, only two gaussians were 
fitted to the histogram in the strip -5 < V (km s~') < 8. The 
parameters of the gaussians yielding the best fit to the data are 
given in Table Q] for the Pleiades strip, Table [2] for the Hyades 
strip, and Table [3] for the Sirius strip. It is worth noting that 
the cr of the background gaussian (33 km/s) is consistent with 
the result listed in table 2 of Famaey et al. (2005), and that the 
parameters of the Hyades and Pleiades gaussians are precisely 
identical (except of course for U and for the absolute number 
of stars present in the strip). 

Table 1. Parameters of the gaussians fitting the U distribution 
of stars from Famaey et al. (2005), in the V-strip corresponding 
to the Pleiades, -28 < V (km s~') < -18. The parameters ki 
give the relative fractions iV;/Mot> where N; is the total number 
of stars in component z, and N tot = 1176 is the total number 
of stars in the V-strip. The fraction f s is the relative fraction of 
stars from the Pleiades moving group with respect to all stars in 
the restricted ±1.4cr u range around the gaussian average, cor- 
responding to the Pleiades zone where the wavelet coefficients 
are non-zero at scale 4 (see Sect. 3). 

background Pleiades 

{/(kms- 1 ) -2 -16 

a u (km s~') 33 5 

k 0.8 0.12 

f s 0.46 



Table 2. Same as Table 1 in the V-strip corresponding to the 
Hyades, -24 < V (km s~') < -12, with iV tot = 1487. 

background Hyades 

U (km s- 1 ) -2 -32 

<Tu (km s~') 33 5 

k 0.8 0.12 

fs 0.52 

The relative contribution of each moving group with re- 
spect to the background at their corresponding location in the 
[U, V] -plane is then obtained by integrating these gaussians 
over a range in U corresponding to ±1.4cr. This has been cho- 
sen to match the regions where the wavelet coefficients are non- 
zero (precisely in the case of the Pleiades, approximately in the 



Table 3. Same as Table 1 in the V-strip corresponding to the 
Sirius moving group, -5 < V (km s~') < 8, with N tot = 1071. 

background Sirius 

f/(kms-') -2 5 

<Tu (km s~') 33 5 

k 0.9 0.1 

L 0.35 



case of the Hyades and Sirius, see Fig. [TJ, i.e. U e [-23, -9] 
for the Pleiades, [-39,-25] for the Hyades, and [-2, 12] for 
Sirius. We find that, in these regions of velocity space, the per- 
centages of stars physically belonging to the Pleiades, Hyades 
and Sirius moving groups are respectively 46% (11 1/239 stars 
in the Pleiades), 52% (150/290 stars in the Hyades) and 35% 
(60/169 in Sirius). 

4.2. Isochrone analysis in parallax space 

We plot in HR diagrams (Fig. 3) the stars kinematically as- 
sociated with the various groups (i.e. stars in the regions de- 
fined by the V-strip and within 1.4cr of the mean of the U- 
gaussian hereabove), together with the isochrone of the ass- 
ociated cluster (100 Myr for the Pleiades, 600 Myr for the 
Hyades, and 300 Myr for Sirius). Let us stress again that these 
HR diagrams mix stars from the background and from the mov- 
ing groups/clusters. Because we perform an analysis in paral- 
lax space hereafter, we chose to plot the (biased) stellar ab- 
solute magnitude My T obtained directly from the inversion of 
the Hipparcos parallax and from the Hipparcos magnitude Hp 
and the Hp - Vj color index, instead of the unbiased abso- 
lute magnitude of Famaey et al. (2005). The absorption cor- 
rection Ay taken from Famaey et al. (2005) has however been 
applied. Note that those absolute magnitudes plotted are thus 
subject to the Lutz-Kelkerbias but that, again, this is not a prob- 
lem since the comparison between isochrone and star locations 
will be made in the parallax space. Finally, the V — Ic index, 
taken as well from Famaey et al. (2005), is derived from the 
Hp - Vt index using the colour transformation from Platais 
et al. (2003). Isochrones from the basic grid (labelled 'c') of 
Schaller et al. (1992) (metallicity Z = 0.02, standard mass 
loss, core overshoot, OPAL opacities) have been used, as re- 
processed by Lejeune & Schaerer (2001) to provide the photo- 
metric indices from the Johnson UBV and Cousins RI bands. 

We then computed for each star in the kinematical samples 
the parallax nT; so expected if the star would lie along the asso- 
ciated cluster isochrone. For the region where RGB and AGB 
overlap, there are more stars in the He-clump or on the RGB, 
since the lifetime is longer in these two phases than on the 
AGB. Therefore, the expected parallax is computed from those 
parts of the isochrones corresponding to the He-clump and 
RGB. Note that, since the isochrones of Lejeune & Schaerer 
(2001) do not cover the whole range of observed V — Ic colors, 
the number of stars to which an "isochrone parallax" could be 
assigned is slightly smaller than the total number of stars in the 
region of velocity space considered. 
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Fig. 2. Top: U distribution of stars from the Famaey et a!. (2005) sample for the Pleiades strip 
-28 < V (km S ) < -18. The thick solid line is a fit to the observed distribution, resulting from the 
superposition of three gaussians (with parameters listed in TableQ, corresponding to the background 
component (thin solid line), to the full Pleiades component and partial Hyades component (thin dashed 
lines). Middle: same as top for the Hyades strip -24 < V (km s _1 ) < -12 (with parameters listed in 
Table [2l. Bottom: same as top and middle for the Sirius strip -5 < V (km s _1 ) < 8 (with parameters 
listed in Table[|] 



Fig. 3- Top: Comparison of the 100 Myr isochrone (Z = 0.02) from Schalleret al. (1992) andLejeune 
& Schaerer (2001) with the stars kinematically associated with the Pleiades cluster/moving group. Note 
that absolute magnitudes have been derived from a simple inversion of the parallax, and are thus bi- 
ased, which is not a problem for the subsequent analysis done in parallax space. Stars whose parallax 
is within Icr of the isochrone parallax are marked with filled circles while others are marked with open 
circles. Middle: Same for the stars kinematically associated with the Hyades cluster/moving group and 
the 600 Myr isochrone. The effect of metallicity is also shown. Bottom: Same for the stars kinematically 
associated with the Sirius cluster/moving group and the 300 Myr isochrone. 
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Fig. 4. Histogram of 6 (see Eq. 7) for the Pleiades moving 
group, together with the gaussian N (see Eq. 8) expected for 
the Pleiades moving group if all stars in the group were on the 
100 Myr isochrone. The subtraction of N from the observed 
histogram is negative (significantly with respect to Poissonian 
errors) at 6 = 0. 

The normalized difference between the "isochrone paral- 
lax" ?ir iso and the measured Hipparcos parallax mmp 



whip - ^is 



(7) 



should then be a random variable following a gaussian of mean 
and standard dispersion 1 for a sample of coeval stars falling 
on the isochrone. Knowing the fraction of stars f s that make 
up the overdensity superimposed on the background in each of 
the regions considered in velocity space (46% in the Pleiades 
zone, 52% in the Hyades zone, and 35% in the Sirius zone), 
we can overplot the corresponding gaussian histogram on the 
total histogram of 6 in the region: this gaussian histogram will 
be equal to 

N = N x A s x f s x MO, 1), (8) 

where N is the total number of stars to which an isochrone 
parallax has been assigned in the region considered and A s is 
the width of the bins of the histogram. 

Then, if a moving group is solely associated with the evap- 
oration of its associated cluster, the subtraction of N from the 
total histogram of 6 in the velocity-space region should yield 
the typical distribution of 6 for the background population. If 
this is not the case, the moving group will be proven to be, at 
least partially, linked with a dynamical (resonant) mechanism. 

4.2.1. The Pleiades moving group 

The histogram of 6 for the 100 Myr isochrone has been plot- 
ted on Fig. 4, together with the gaussian histogram N corre- 
sponding to the hypothesis that all stars of the Pleiades mov- 
ing group (i.e. 46% of the stars in the Pleiades velocity-zone) 
are falling on the 100 Myr isochrone. Very clearly, the actual 
amount of stars really present in the ^-interval covered by the 



Fig. 5. Same as Fig. 4 for the Hyades moving group and the 
600 Myr isochrone. 



gaussian is much too low as compared to N to be compatible 
with the hypothesis that 46% of the stars are coeval at 100 Myr. 
If one would subtract N to the actual histogram to find the back- 
ground distribution in parallax space, one would find that the 
background density distribution is highly negative close to the 
100 Myr isochrone, which would of course not make any sense. 

This is thus a robust proof that the Pleiades moving group 
is not associated with the evaporation of the Pleiades cluster. 
Note that the excess of stars in the 6 histogram compared to 
the expectation for an evaporated cluster is more prominent for 
positive values of 6: these stars are closer (and thus intrinsi- 
cally fainter) than what they should be to fall on the 100 Myr 
isochrone, meaning that correcting for the Malmquist bias (i.e. 
the fact that intrinsically faint objects are missing from such 
a magnitude-limited sample) would make the problem even 
worse. Actually, a qualitative look at Fig. 3 reveals that 15 stars 
are clearly falling close to the 100 Myr isochrone, and are thus 
likely members of the evaporated Pleiades cluster. 

4.2.2. The Hyades moving group 

The histogram of 6 for the 600 Myr isochrone has been plotted 
on Fig. 5, together with the gaussian histogram N correspond- 
ing to the hypothesis that all stars of the Hyades moving group 
(i.e. 52% of the stars in the Hyades velocity-zone) are falling 
on the 600 Myr isochrone. Very clearly, subtracting N from the 
observed histogram would imply that the background popula- 
tion on which the Hyades group is superimposed has a gap in 
parallax space close to the 600 Myr isochrone. 

To check whether this is the case, we computed the 6- 
histogram in a comparison box symmetrical to the Hyades with 
respect to U = -2, i.e. U = [21,35], V = [-24,-12]. This 
histogram (grey-shaded in Fig. 6) should be representative of 
the ^-distribution of the background population in the Hyades 
velocity-zone. We see in Fig. 6 that it is in disagreement with 
the ^-distribution in the Hyades velocity-zone after subtract- 
ing N (Fig. 6 bottom; normalized to the number of stars in the 
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Fig. 6. The shaded area is the histogram of 6 for the 600 Myr 
isochrone in the comparison box U = [21,35], V = [-24,-12], 
compared with (top) the histogram of 6 in the Hyades box 
(background+moving group), and (bottom) the histogram of 5 
expected for the background if all stars of the Hyades overden- 
sity were coeval (600 Myr), after subtracting the N histogram 
of Fig. 5. All histograms are normalized to the number of stars 
in the comparison box. 



comparison box), but compatible with the full ^-distribution in 
the Hyades velocity zone (Fig. 6 top). We tested the influence 
on our result of the metallicity adopted for the isochrone (see 
Fig. 3), and of the actual age assigned to the Hyades cluster (by 
trying 800 Myr instead of 600), and found the same disagree- 
ment. We thus conclude that it is impossible that the Hyades 
moving group, representing 52% of the stars in the Hyades 
velocity-zone, is entirely made of stars evaporated from the 
Hyades cluster. 

4.2.3. The Sirius moving group 

The histogram of 5 for the 300 Myr isochrone has been plotted 
on Fig. 7, together with the gaussian histogram N correspond- 
ing to the hypothesis that all stars of the Sirius moving group 
(i.e. 35% of the stars in the Sirius velocity-zone) are falling on 
the 300 Myr isochrone. 

We then computed (grey-shaded in Fig. 8) the 5-histogram 
in a comparison box symmetrical to Sirius with respect to 
U = -2, i.e. U = [-16,-2], V = [-5,8]. This histogram 
should be representative of the ^-distribution of the background 
population in the Sirius velocity-zone. We compared it to both 
the full ^-distribution in the Sirius velocity zone, and to the 5- 
distribution in the Sirius velocity-zone after subtracting N. In 
each case the histogram was normalized to the number of stars 
in the comparison box. Once again, while the observed his- 
togram of S in the comparison box is compatible with the to- 
tal one in the Sirius box, it is in Poissonian disagreement with 
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Fig. 7. Same as Figs. 4 and 5 for the Sirius moving group and 
the 300 Myr isochrone. 
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Fig. 8. The shaded area is the histogram of 6 for the 300 Myr 
isochrone in the comparison box U = [-16,2], V = [-5,8], 
compared with (top) the histogram of 5 in the Sirius box (back- 
ground+moving group), and (bottom) the histogram of 6 ex- 
pected for the background if all stars of the Sirius overdensity 
were coeval (300 Myr), after subtracting the N histogram of 
Fig. 7. All histograms are normalized to the number of stars in 
the comparison box. 



the histogram obtained after subtracting N (especially around 
5-2). However, the disagreement is less obvious than for the 
Hyades: this is because the Sirius overdensity superimposed 
on the background is much less prominent than the Hyades one 
(see Fig. 2). 
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5. Conclusion and perspectives 

We have reanalyzed the data of Famaey et al. (2005) by using 
a wavelet transform in order to identify the main low-velocity 
moving groups on scales smaller than 30 km/s, and by subse- 
quently comparing the location of moving group stars in paral- 
lax space to isochrones of the embedded clusters. 

This leads to the fascinating result that moving groups do 
appear in pairs in the [U, V] -diagram (The Hyades and Pleiades 
groups are similar to the two distinct structures in the Hercules 
stream), which could provide new interesting constraints on 
dynamical models designed to reproduce these features (see 
e.g. figure 9 of De Simone et al. 2004, figure 6 of Quillen & 
Minchev 2005, and figure 10 of Chakrabarty 2007) if their ori- 
gin is in resonant trapping by the bar, spiral arms, or a mix of 
both. It would also be of high interest to check whether some 
differences in metallicities could be observed between the two 
distinct structures of the Hercules stream to check whether they 
could have an origin at different galactocentric radii. 

Moreover, this study confirms the dynamical (resonant) na- 
ture of the Pleiades, Hyades and Sirius moving groups (see 
also Famaey et al. 2007), because the fraction of stars mak- 
ing up each velocity-space overdensity superimposed on the 
background is higher than the fraction of stars compatible 
with the isochrone of the embedded cluster. The similarity be- 
tween the parallax distribution with respect to the isochrone 
in the velocity-zones representative of the background and of 
the moving groups is in fact striking (see Figs. 6 and 8), and 
again demonstrates that the group stars share the same distri- 
bution as the background stars in the HR diagram. These low- 
velocity moving groups have thus a different nature than, e.g., 
the HR1614 group (U * [0, 20], V [-65, -45], see De Silva 
et al. 2007) which was recently proven to be an evaporated clus- 
ter. 

The next check of this resonant origin of the Pleiades, 
Hyades and Sirius moving groups should be a detailed chem- 
ical tagging of stars kinematically associated with the groups 
(Freeman & Bland-Hawthorn 2002) from their element abun- 
dance patterns. For instance, we have identified here 290 giant 
stars kinematically associated with the Hyades, out of which 
150 are making up the Hyades overdensity in velocity space 
(or Hyades moving group). In the resonant scenario, we expect 
less than 150 giant stars in this sample to have chemical abun- 
dances in precise accordance with those of the Hyades cluster 
(De Silva et al. 2006), although we still expect some of them 
to be part of the evaporating cluster. Note that we would how- 
ever expect those 150 stars to be more metal-rich on average 
than the background because of the likely origin of the moving 
group in the inner parts of the Galaxy (see e.g. Daflon & Cunha 
2004 for the galactic metallicity gradient). In the case of the 
Sirius moving group, since the overdensity superimposed on 
the background is much less prominent than the Hyades one 
(see Fig. 2), and since the evaporated Ursa Major cluster could 
still make a large part of the moving group, the real dynamical 
overdensity should not represent as large a part of the sample 
as for the Hyades in such a chemical tagging analysis. Finally, 
another independent check of the resonant origin of the moving 
groups would be a detailed study of the individual kinematics 



of open clusters in the galactic disk (e.g., Frinchaboy 2006): 
within the resonant scenario, open clusters with a wide range of 
ages should show the same overdensities in the [U, V] -diagram 
as the individual stars studied in this paper. 

Acknowledgements. We thank Pavel Kroupa for useful discussions, 
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